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■ ABSTRACT 

We examine the power source of 41 local Ultraluminous Infrared Galaxies (ULIRGs) 
using archival infrared (IR) and optical photometry. We fit the observed Spectral 
Energy Distributions (SEDs) with starburst and AGN components; each component 
being drawn from a family of templates. We find all of the sample require a starburst, 
whereas only half require an AGN. In 90% of the sample the starburst provides over 
i-G ■ half the IR emission, with a mean fractional luminosity of 82%. When combined with 

other galaxy samples we find that starburst and AGN luminosities correlate over 6 
q ■ decades in IR luminosity suggesting that a common factor governs both luminosities, 

?-h ' plausibly the gas masses in the nuclear regions. We find no trend for increasing frac- 

, tional AGN luminosity with increasing total luminosity, contrary to previous claims. 

■ We find that the mid-IR -£7.7/67.7 line-continuum ratio is no indication of the star- 
burst luminosity, or the fractional AGN luminosity, and therefore that .FY.7/C7.7 is not 
a reliable diagnostic of the power source in ULIRGs. The radio flux correlates with the 
starburst luminosity, but shows no correlation with the AGN luminosity, in line with 
previous results. We propose that the scatter in this correlation is due to a skewed 
starburst IMF and/or relic relativistic electrons from a previous starburst, rather than 
contamination from an obscured AGN. We show that most ULIRGs undergo multiple 
starbursts during their lifetime, and by inference that mergers between more than 
two galaxies may be common amo ngst ULIRGs . Our results support the evolution- 
ary model for ULIRGs proposed bv lFarrah et all l)2001f) . where they can follow many 
different evolutionary paths of starburst and AGN activity in transforming merging 
spiral galaxies into elliptical galaxies, but that most do not go through an optical QSO 
phase. The lower level of AGN activity in our local sample than in z ~ 1 HLIRGs im- 
plies that the two samples are distinct populations. We postulate that different galaxy 
formation processes at high-z are responsible for this difference. 

Key words: infrared: galaxies - galaxies: active - galaxies: Seyfert - galaxies: star- 
burst - Quasars: general 
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1 INTRODUCTION 

One of the most important results from the Infrared As- 
tronomical Satellite (IRAS) all sky surveys was the detec- 
tion of a new class of galaxy where mos t of the bolomet- 
ric emission lies in the in frared waveband (ISoifer et a1ll984l: 
ISanders fc Mirabeflll996h . These 'Luminous Infrared Galax- 
ies' (LIRGs), become the dominant extragalactic population 



at IR luminosities above 10 Lq, with a higher space density 
than all other classes of galaxy of comparable bolometric lu- 
minosity. At the most luminous end of the IRAS galaxy pop- 
ulation lie the Ultraluminous Infrared Galaxies (ULIRGs), 
those with Lt r > 1O 12 L0. Although ULIRGs are rare in 
the local Universe, the ir luminosity function shows st rong 
evolution with redshift jVgjlleirx^Sandi^ and 
deep sub-mm surveys llBarger^^i]F99^]H^gh^set^ill993 : 
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Table 1. All Ultraluminous Infrared Galaxies at z < 0.1 



IRAS Name 


Other names 


RA (2000) 
hh:mm:ss 


Dec 

o / // 


z 


Optical spectrum 




00198-7926 




00 21 52.9 


-79 10 08.0 


0.0728 


Sy2 


0.370 


00199-7426 




00 22 07.0 


-74 09 41.6 


0.0964 


- 


0.078 


00262+4251 




00 28 54.2 


+43 08 15.3 


0.0927 


LINER 


0.112 


00335-2732 




00 36 00.5 


-27 15 34.5 


0.0693 


Starburst 


0.147 


01388-4618 




01 40 55.9 


-46 02 53.3 


0.0903 


HII 


0.121 


02364-4751 




02 38 13.1 


-47 38 10.5 


0.0983 


- 


0.063 


03068-5346 




03 08 21.3 


-53 35 12.0 


0.0745 


- 


0.057 


04232+1436 




04 26 05.0 


+14 43 38.0 


0.0796 


LINER 


0.113 


05189-2524 




05 21 01.5 


-25 21 45.4 


0.0426 


Sy2 


0.252 


06035-7102 




06 02 54.0 


-71 03 10.3 


0.0795 


HII 


0.112 


06206-6315 




06 21 01.2 


-63 17 23.2 


0.0924 


Sy2 


0.074 


08572+3915 




09 00 25.4 


+39 03 54.4 


0.0584 


LINER 


0.229 


09111-1007 




09 13 38.8 


-10 19 20.3 


0.0541 


HII+Sy2 


0.067 


09320+6134 


UGC 05101 


09 35 51.7 


+61 21 11.3 


0.0394 


LINER+Sy2 


0.090 


09583+4714 




10 01 31.2 


+46 59 44.0 


0.0859 


Syl+Sy2 


0.182 


10035+4852 




10 06 46.1 


+48 37 44.1 


0.0648 


Starburst 


0.062 


10190+1322 




10 21 42.5 


+13 06 53.9 


0.0766 


HII 


0.114 


10494+4424 




10 52 23.6 


+44 08 47.3 


0.0921 


LINER 


0.047 


10565+2448 




10 59 18.1 


+24 32 34.3 


0.0431 


HII+LINER 


0.094 


12112+0305 




12 13 45.7 


+02 48 40.3 


0.0733 


LINER 


0.060 


12540+5708 


Mrk 231 


12 56 14.2 


+56 52 25.2 


0.0422 


Syl 


0.271 


13428+5608 


Mrk 273 


13 44 42.1 


+55 53 12.6 


0.0378 


Sy2 


0.105 


14348-1447 




14 37 38.4 


-15 00 22.8 


0.0827 


LINER 


0.072 


14378-3651 




14 40 58.9 


-37 04 33.0 


0.0676 


LINER+Sy2 


0.085 


15250+3609 




15 26 59.4 


+35 58 37.5 


0.0552 


LINER 


0.182 


15327+2340 


Arp220 


15 34 57.1 


+23 30 11.5 


0.0181 


LINER 


0.076 


17132+5313 




17 14 20.5 


+53 10 30.4 


0.0509 


HII+AGN 


0.109 


17208-0014 




17 23 21.9 


-00 17 00.4 


0.0428 


HII 


0.053 


18470+3233 




18 48 54.2 


+32 37 31.0 


0.0784 


HII 


0.104 


19254-7245 


Super-Antena 


19 31 21.6 


-72 39 21.7 


0.0617 


Sy2 


0.226 


19297-0406 




19 32 20.7 


-04 00 06.0 


0.0857 


HII 


0.084 


19458+0944 




19 48 15.7 


+09 52 05.0 


0.0999 




0.066 


20046-0623 




20 07 19.4 


-06 14 26.0 


0.0844 


Starburst 


0.141 


20414-1651 




20 44 18.2 


-16 40 16.2 


0.0871 


HII 


0.079 


20551-4250 




20 58 26.9 


-42 39 06.2 


0.0428 


HII 


0.149 


21130-4446 




21 16 18.5 


-44 33 37.7 


0.0926 


HII 


0.048 


21504-0628 




21 53 05.5 


-06 14 49.9 


0.0776 




0.111 


22491-1808 




22 51 49.3 


-17 52 23.4 


0.0778 


HII 


0.101 


23128-5919 




23 15 47.0 


-59 03 16.9 


0.0446 


HII+Sy2 


0.147 


23365+3604 




23 39 01.3 


+36 21 08.7 


0.0645 


HII+LINER 


0.114 


23389-6139 




23 41 43.6 


-61 22 50.9 


0.0928 


HII 


0.067 



Coordinates, optical spectral classifications and redshifts were taken from the NASA Extragalactic Database. "Ratio of the IRAS 
25^m flux to the IRAS 60/im flux. 



Conversely, many ULIRGs display em ission lines character- 
istic of Seyferts llSanders et alHl988T). An early evolu tion- 
ary model for ULIRGs was that of lSanders et all lll98ct) who 
suggested, based on the similar space density, bolometric 
emission and luminosity function of ULIRGs and QSOs in 
the local Universe, that ULIRGs are the dust enshrouded 
precursors to optical QSOs and that all QSOs emerge from 
a luminous infra red phase. Conver sely, a more recent evo- 
lutionary model dFarrah et alll200lT) proposes that ULIRGs 
are not a simple dust shrouded precursor to optical QSOs 
but instead follow multiple evolutionary paths. 

In this paper we examine the power source in ULIRGs, 
their evolution, and their relationship to high-z IR luminous 
galaxies using archival photometry for a sample of 41 local 
ULIRGs, and advanced radiative transfer models for star- 
bursts and AGN. We also examine the origin of the radio- 
IR correlation in ULIRGs and the power of mid-IR spec- 



lEales et ai200rilScott et a^OOi&rFax et al2002h have found 
that systems with ULIRG-like luminosities are very nu- 
merous at z ^ 1. ULIRGs are thus an important popula- 
tion in understanding the cosmic history of star formation. 
ULIRGs are also thought to play a role in the evolution 
of spiral and elliptical galaxies. Nearly all ULIRGs are ob- 
served to be ongoing mergers between two or more spirals 
JSanders et al Il988l : iFarrah et ajl200ll) . and it is thou ght 
that such a merger will form an elliptical dBarnesHl989lL 

Despite this, the evolution of ULIRGs, their power 
source and the trigger behind the IR emission are 
poorly understood. Although it is now accepted that 
a mixture of star formation and AGN activity powers 
the IR emission, the dominant power source, and how 
ULIRGs evolve, are unknown. Ther e are similarities between 
ULIRGs and starburst galaxies (iJoseph fc Wright Il985t 
iRowan-Robinson fc Crawford! Il989l : ICondon et all Il99l¥ . 
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troscopy as a diagnostic of the active power source. Sample 
selection and data analysis are described in *J5|and Re- 
sults are presented in Jl]and notes on individual objects are 
given in ^| Discussion is presented in SjfJ and conclusions 
are summarized in J7| We adopt Ho = 65 km s _ Mpc - , 
f2o = 1-0 and A = 0.0 and quote all luminosities in this 
system. Unless otherwise stated, the terms 'IR luminosity', 
'starburst luminosity', and 'AGN luminosity' refer to the 
luminosity over the wavelength range 1 — 1000/im. Lumi- 
nosities are quoted in units of bolometric solar luminosities, 
where L B = 3.826 x 10 26 Watts. 



2 THE SAMPLE 

We assembled from the literature a sample of objects com- 
prising all ULIRGs with z ^ 0.1. This was complicated by 
the fact that different aut hors use different m inimum lu- 
minosities for_JJL2RGs_^gJCorido^ L40-120 =S 
10 1102 L^. lciements et all jl996h : L«n ~0 11 - 77 L S ). We first 
adopted our own minimum IR luminosity for a ULIRG, 
namely a rest-frame 1 — 1000/im luminosity (hereafter re- 
ferred to as Li r ) of greater than 10 12 L©. Using a luminosity 
derived across a broad wavelength range is to be preferred, as 
luminosities derived across smaller wavelength ranges may 
erroneously exclude or include objects due to their contin- 
uum shape. We then assembled a parent sample of all ob- 
jects claimed to be ULIRGs, or close to ULIRG luminosity, 
at z 0.1, from the IRAS Point Source Catalogue, the 
IRAS Fai nt Source C atalogue, and the sample presented by 
IClements et all lll996l) . We then determined their 1 — 1000/im 
luminosities using the methods described in sections [3] and 
2] and excluded those objects that did not satisfy our lu- 
minosity criterion. From a parent sample of 73 objects, 32 
were excluded. The final sample of 41 objects, presented in 
TableQ are all contained within the IRAS Point Source Cat- 
alogue. This final sample thus comprises a complete sample 
of ULIRGs at z < 0.1 with Li r > 1O 12 L0, and is thus free 
from selection effects arising from infrared colours, or lumi- 
nosities measured across smaller wavelength ranges. 

Fluxes were assembled from online catalogues and 
from the literature. Optical and near-IR fluxes were 
taken from the APM and 2MASS datab a ses a nd also 
from ICarico et all Jl98Sl) & ISpinoglio et"al lll995D . IRAS 
fluxes were taken from the IRAS Faint Source Survey. 
The XSCANPI software was used to derive IRAS fluxes 
where only upper limits were presen t in the cat a logues . 
Other infrared data were taken from ICarico et all il988}): 
iMaiolino et all dl995l): iKlaas et al <1997lk iRigopoulou et a 
|l999[ l: iDale et al feOOOl): IZink et all feOOCft ; iTran et a 
fcOOlli fc IKlaas et all (120011). Sub-millimetre data were 
taken from IRigop oulou. La wrence fc Rowan- Robin son! 
l|l996ft : iDunne et al fcOOOl) fc iDunne fc Ealesl feOOll) . The 



3 INFRARED EMISSION MODELS 



3.1 Starburst Models 



assembled optical and near-infrared photometry is listed in 
Table HJ. 



To model the IR 
activity we used 



emission due 
the starburst 



to starburst 
models of 

lEfstathiou. Rowan-Robinson fc Siebenmorgenl i200(J) . 
These models consider an ensemble of evolving HII re- 
gions containing hot young stars, embedded within Giant 
Molecular Clouds (GMCs) of gas and dust. The com- 
po sition of the dust is giv e n by the dust grain model 
of ISiebenmorgen fc KruegeH (Il992f) . which includes the 
Polycyclic Aromatic Hydrocarbons (PAHs) thought to be 
responsible for the 7.7/^m emission feature. The stellar 
populations within the GMCs evolve according t o the 
stellar synthesis codes of iBruzual fc Charlotl 1^993) . The 
star formation rate is assumed to decay exponentially with 
an e-folding timescale of 2 x 10 7 years. The models vary 
in starburst age from zero years up to 7.2 x 10 7 years (11 
discrete values) and in effective visual optical depth from 
tv — 50 to tv = 200 (4 discrete values). A plot showing 
the range in starburst SEP template shapes can be found 
in figu re 3 of Efstat hiou. Rowan-Robinson fc Siebenmorgenl 
(2000), with additional information in figure 1 of the same 
paper. 



3.2 AGN Torus Models 

To model t he IR emission due to AGN, we used t he AGN 
models of lEfstathiou fc Rowan- Robinson! jl995l) . These 
models incorporate accurate solutions to the axially sym- 
metric radiative-transfer problem in dust clouds to model 
the IR emission from dust in active galactic nuclei. Dust 
composition is g i ven b y the multigrain dust model of 
iRowan-Robinsonl lll992l) . We have used the thick tapered 



1 A table containing the assembled infrared 
sub-millimetre photometry can be obtained 
http: / /spider. ipac.caltech.edu/staff/duncan/irphotom.tex 



and 
from 



disk models following an r _1 density distribution, as this 
subset of models has been found to be most successful in 
satisfying the observational constraints of AGN. The AGN 
models vary in torus opening angle from 0° to 90° (15 
discrete values) and in UV equatorial optical depth, from 
tuv = 1000 to tuv = 1500 (3 discrete values). A plot show- 
ing the range in A GN torus SEP template shape s can be 
found in figure 5 of lEfstathiou fc Rowan-Robinsonl (^995). 



4 RESULTS 

4.1 Spectral Energy Distributions 

We combined the archival photometry from the optical to 
the millimetre to fit Spectral Energy Distributions (SEDs) 
for each object. Goodness of fit was examined by using the 
reduced x 2 statistic. Fits for all sources were good, with 
XfLst ^ 5 in all cases, and xLst ^ 4 in all but one case 
(this object, IRAS 17208-0014, is discussed further in Sec- 
tion |2J • To determine the errors on the model parameters 
we explored reduced x 2 space between xLst an d Xbest + 2. 
The resulting luminosities, model parameters and errors are 
given in Table [2] The SEDs are presented in Figure 

Emission from unobscured population II stars from the 
merger progenitors is not included in either the starburst 
or AGN models, and care must be taken in accounting for 
this, as a recent Hubble Space Telescope (HST) study of 
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Table 2. Ultraluminous Infrared Galaxies: Optical and near-IR data 



Name 


U 


B 


V 


R 


J 


H 


K 


L 


00198-7926 










619 + 41 


764 ± 51 


565 ± 38 


232 ± 15 


00199-7426 










1.0 + 0.1 


5.4 + 0.6 


6.4 + 0.7 




00262+4251 










5.4 + 1.6 




5.5 ± 1.6 




00335-2732 








0.94 + 0.1 


4.2 + 0.4 


3.3 ± 0.3 


5.0 + 0.5 




01388-4618 




0.31 ± 0.04 














02364-4751 




0.44 + 0.05 














03068-5346 


















04232+1436 










5.5 + 1.6 




10.6 + 3.2 




05189-2524 




0.66 ± 0.07 


0.65 + 0.07 


0.76 ± 0.08 


14.7 + 0.9 


30.8 ± 2.1 


57.1 + 3.8 


136 ± 30 


06035-7102 










2.1 + 0.2 


2.8 ± 0.3 


4.2 + 0.4 




06206-6315 










1.3 + 0.1 


1.7 + 0.2 


2.8 + 0.3 




08572+3915 




0.51 ± 0.05 


0.65 + 0.07 


0.90 + 0.1 


1.7 + 0.1 


3.0 + 0.2 


3.9 + 0.3 


50 + 10 


09111-1007 










3.4 + 0.3 


5.5 ± 0.6 


7.4 + 0.7 




09320+6134 




0.96 + 0.1 


10.91 ± 1.2 


1.11 + 0.1 


12.9 ± 0.9 


21.2 ± 1.4 


29.2 + 1.9 


33 + 7 


09583+4714 










1.2 + 0.1 


1.5 ± 0.2 


2.0 + 0.2 




10035+4852 


















10190+1322 










2.3 + 0.2 


3.0 + 0.3 


2.1 + 0.2 




10494+4424 










1.0 + 0.1 


2.0 + 0.2 


2.8 + 0.3 




10565+2448 










14.0 + 1.4 


17.0 + 1.7 


18.0 + 1.8 




12112+0305 




0.16 ± 0.02 


0.25 + 0.03 


0.32 + 0.03 


2.0 + 0.2 


2.9 ± 0.3 


3.5 + 0.4 


2.2 + 0.5 


12540+5708 


2.77 + 0.5 


7.34 + 1.24 


13.1 ± 2.2 


11.5 + 0.1 


49 + 3 


107 + 7 


186 ± 12 


368 ± 72 


13428+5608 




0.74 + 0.07 


1.03 ± 0.1 


1.76 + 0.2 


9.0 + 0.9 


13.5 ± 1.4 


16.1 ± 1.6 


17.4 + 1.7 


14348-1447 




0.22 ± 0.02 


0.24 + 0.02 


0.43 ± 0.04 


1.7 + 0.2 


2.5 ± 0.3 


3.3 + 0.3 




14378-3651 










2.5 + 0.3 


2.8 ± 0.3 


5.0 + 0.5 




15250+3609 




0.26 ± 0.03 


0.27 + 0.03 


0.42 ± 0.04 


3.1 + 0.3 


4.0 + 0.4 


4.1 + 0.4 


4.0 + 0.4 


15327+2340 




0.96 + 0.1 


3.16 ± 0.3 


4.07 + 0.4 


9.6 + 1.0 


17.9 ± 1.8 


23.0 + 2.3 


22.7 + 2.3 


17132+5313 










5.6 + 0.6 


8.2 ± 0.8 


8.9 + 0.9 


10.5 + 1.6 


17208-0014 


















18470+3233 










2.7 + 0.4 




3.6 + 0.5 




19254-7245 


















19297-0406 


















19458+0944 


















20046-0623 










3.9 + 0.6 




4.0 + 0.6 




20414-1651 


















20551-4250 




3.93 ± 0.34 




5.04 + 0.34 


22.9 ± 5.9 


30.2 ± 7.8 


25.7 + 6.7 




21130-4446 








0.89 + 0.09 










21504-0628 










5.5 + 0.8 




5.7 + 0.9 




22491-1808 




0.22 ± 0.02 


0.20 + 0.02 


0.26 + 0.03 


2.1 + 0.2 


2.7 + 0.3 


2.6 + 0.3 


2.4 + 0.2 


23128-5919 




4.09 + 0.4 




7.69 + 0.7 


9.6 + 0.6 


12.6 + 0.8 


11.9 + .08 




23365+3604 










7.5 + 0.4 




6.5 + 0.3 





23389-6139 0.35 ± 0.04 



All fluxes are given in mjy and are in the observed frame, la errors are quoted. 



ULIRGs jFarrah et, afeOOlf ) found that the optical emission 
was in most cases dominated by old stellar populations and 
not by light from a starburst or AGN. The range in possible 
optical/near-IR SEDs for evolved stellar populations is very 
large, and in addition an obscured starburst or AGN can in 
principle contribute significantly in the optical or near-IR. 
Constraining the properties of the evolved stellar compo- 
nent in any of our sample with the limited available optical 
photometry therefore proved impossible. We have therefore 
assumed that all emission shortwards of 3.5/im in the rest 
frame of the objects contains a significant but unquantified 
contribution from old stellar populations. As this contribu- 
tion could lie between 0% and 100%, any measured flux at 
a rest-frame wavelength shorter than 3.5/im is treated as a 
3a upper limit in the fitting. 



4.2 Star Formation Rates 

Estimating obscured star formation rates from IR data is 
based on silicate and graphite dust grains absorbing the 
optical and UV light from young stars and reradiating in 
the IR and sub-mm. A recent estimate for deriving star 
formation rates from IR lu minosities has been made by 
Rowan- Robinson et a 

M,. aii =2.6xl0- 10 ^ (1) 

where M*, a ii is the star formation rate, e is the fraction of 
optical/UV light from the starburst that is absorbed by dust 
and re-emitted in the IR, and Lm is the 60/^m starburst lu- 
minosity. The factor <j> incorporates the correction between 
a Salpeter IMF and the true IMF (xl.O for a Salpeter 
IMF, x3.3 for a Miller-Scalo IMF), and a correction for 
the assumed upper and lower stellar mass bounds for the 
stars forming in the starburst (xl.O if the mass range is 
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0.1 < Mq < 100, x 0.323 if the mass range spans only OBA 
type stars, i.e. 1.6 < Mq < 100). We have assumed that 
all of the optical/UV light from the starburst is absorbed 
(e ~ 1.0), and that the IMF of the starburst is a Salpeter 
IMF forming stars across the mass range 0.1 < Mq < 100, 
(<j> = 1.0). The star formation rates for our sample calculated 
using Equation are given in Table[31 



4.3 Dust Masses and Temperatures 

The starburst and AGN models do not assume a monolithic 
dust temperature, but instead invoke a range of dust tem- 
peratures spanning 10K to 1000K. An estimate of the mass- 
weighted mean dust temperature can be obtained by fitting 
an optically thin greybody function of the form: 



S v = v f5 B{v,T dus 



(2) 



to the FIR SED over the wavelength range 200 — 1000/^m, 
which measures the colour temperature of the system. In 
Equation S v is the source flux at a frequency v, T dust 
is the dust temperature, B(y, Tdust) is the Planck function, 
and P is the frequency dependence of the grain emissivity. 
Monolithic dust temperatures derived in this way are un- 
physical simplifications in AGN, but can serve as a useful 
comparison with previous work. Mass-weighted mean dust 
temperatures for our sample are presented in Table |H] 

Gas and dust masses can be estimated in two ways. The 
first way is to calculate the gas mass directly from the mod- 
els using equation 4 of iFarrah et all (|2002b). This however 
gives the total gas mass at the current age of the starburst 
without accounting for gas that is blown out of the nuclear 
regions by supernovae or superwinds, and thus may over- 
estimate the amount of gas and dust in the system. The 
second approach is to estimate the dust mass directly from 
the sub-mm SED. An approach for estimating dust and gas 
masses in galaxies using sub-mm data is to assume the sys- 
tem is opticall y thin at the s e wav elengths, and to follow the 
prescription of iHildebrandl lll983l) : 



M d 



1 + 2 K(v r )B(l> r , Tdust) 



(3) 



where v and v r are the observed and rest frame frequen- 
cies respectively, S Vo is the flux in the observed frame, 
B(u r , Tdust) is the Planck function in the rest frame and 
Tdust is the dust temperature. The gas mass is then ob- 
tained by assuming a fixed gas to dust ratio. For the most 
extreme IRAS galaxies, th e best current estimate of the ga s 
to dust ratio is 540 ± 290 JSanders. Scoville fc Soiferjll99ll) . 
For comparison, the gas to dust ratio in spira l galaxies is 
thought to b e ~ 500 iPevereux fc Ypundll99(il) . and ~ 700 
in ellipticals JWiklind fc Henkeil995lT The mass absorption 
coefficient in the rest frame, n r , is taken to be: 



K r = 0.067 



( * V 

V2.5 x 10 11 / 



(4) 



in units of m 2 kg 1 . Dust masses and temperatures calcu- 
lated using Equations |21 and |3 are given in Table [3] 



5 INDIVIDUAL OBJECTS 

In this section we present further discussion on objects in 
our sample that are either well-studied by previous authors, 
or that show interesting or unusual features. 

00198-7926: This source is a double nucleus system with a 
large tail, with a Sy2 optical spectrum. The second nucleus 
is probably a large region of extranuclear unobscured star 
formation, an d there are a number of HII regions embedded 
in the tail jHeisler fc VadeJ Il994f) . A deep Bepposax 
observation gave no detection, giving an uppe r limit to the 
X-ray flux of F 2 -io < 1 X 10~ 13 erg cm -2 s" 1 iRisaliti et all 
2000). From our results, the properties of this object mark 
it as unusual. We find that the IR emission from this object 
is mostly starburst in origin. The starburst is comparatively 
old, and can thus account for both the sub-mm emission 
and all of the unobs cured star formation observed by 
iHeisler fc Vadeil Jl994f) . We interpret the 'warm' infrared 
colour of this object as being due to the old starburst 
rather than an AGN. There is however comparatively little 
IR data available for this object and the resulting AGN 
upper limit is weak. The best fit model, shown in Figure 
includes an AGN component. 

05189-2524'- This object appears to be a late-stage merger. 

It possesses a single, compact, v ery red nucleus with a 

Sy2 spectrum JVeillcux ct al 1995), bisect ed by a dust lane 

l|Scoville et al boOOl) . lYoung et all dl996Tl observed broad 

lines in polarized flux in this object, suggesti ng the presence 

of an obscured AGN. Later observations JPudlevI 119991) 

observed the 11.3/xm dust feature, suggesting that this 

object also contains a buried starburst. The X-ray spectrum 

( Risal itT et all [2000) is well fitted by a two-component 

1 »q+ - 35 
1 -° 3 -o .34 

and 

a thermal component with kT = 0.88jlo'|§ keV. Overall the 
power law component can be interpreted as arising from a 
Compton thin AGN, with the thermal component either 
due to the AGN or to starburst activity. For our SED fitting 
we compiled an a dditional N band (10. 6/jm) flux of 498mJy 
± lOOmJy from iMaiolino et all ((1995) . From our results, 
we find that this object contains both a starburst and an 
AGN, where the AGN provides a significant fraction of the 
total IR luminosity. 

09320+6134: This object posses ses a single bright nucleus 
and a disturbed spiral structure llScoville et al 2000). There 
is one large tail extending ~ 38kpc, and a second tail that 
forms almost a complete ring with a total length of ~ 65kpc. 
The morphology has been interpr eted as a late stage merger 
between two large spiral galaxies llSanders et al|l 988|). possi- 
bly in volving a third gas-poor dwarf galaxy jMaiewski et all 
1993). The optical spectrum is a combination of a LINER 
and a Sy2 jGoncalves, Veron-Cetty J _fc_V^ . Near- 

IR spectroscopy by Imanis hi. Dudley fc Malonevl (l200lT) 
shows that this obj ect is likely to contain a h e avily obscured 
AGN. Conversely, iLutz. Veilleux fc Genzel (I1999T) classify 
this object as starburst powered on the basis of mid-IR spec- 
troscopy, and th e sub-mm emission is also consistent with 

a starburst iRigopoulou. Lawrence fc Rowan- Robinsonl 

1996). We compiled an additional 350/xm upper limit of 
2664m Jy, and an additional 800/jm flux of 143mJy ±25mJy 
from iRigopoulou. Lawrence fc Rowan-Robinsonl Jl996l) . 
We find that this system is a composite object containing 



model, consisting of a power law with F — _l.o»_ 34 , 
absorbed by a column density of Nb ~ 4.7ti't cm -2 , 
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Figure 1. Best fit Spectral Energy Distributions for the 41 ULIRGs in our sample. In each case the solid line is the combined best- fit 
model, the dotted line is the Starburst component and the long dashed line is the AGN component. 



both a starburst and an AGN, a result that is consistent 
with previous observations. From the SED fit presented in 
Figure it can be seen that the AGN has a higher near-IR 
flux than the starburst component, but that the AGN 
contribution in the mid-IR is small, and is negligible in the 
sub-mm. This object is a case in point that, in order to 
understand the power source in ULIRGs it is necessary to 
have data spanning a wide wavelength range. 



10565+2448: We compiled an additional 350/xm flux of 
1240m Jy ±~248mJy, and a 750/mi flux of 85mJy ±17mJy 
from iRjgoDoulou. Lawrence fc Rowan-Robinsonl l)l99fl) . 
The AGN in this object, although weak compared to the 
starburst, is required for an acceptable SED fit. 
12112+0305: This syste m contains two se parate nuclei 
with a pair of tidal tails llScoville et a]ll200nh. The optical 
spectrum is that of a LINER llVeilleux. Sanders fc Kirra 
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Figure 1 - continued 



1999). This source is classi fied as a starburst on the 
basis of its mid-IR spectrum iRjgonoulou e t al llflflflh . We 
compiled an additional 800/xm flux of 50mJ y ±13mJy from 
iRigoDoulou. Lawrence fc Rowan-Robinsonl dl996f) . The 
extensive IR data available for this object allow us to place 
strong constraints on the power source. We find that this 
object is powered by a starburst, with a severe upper limit 
on any AGN contribution. 

12540+5708: This galaxy is commonly known as Mrk 231, 



and has a Syl optical spectrum. The IR luminosity is 
thoug ht to be powered by both s tarburst and AGN a ctivity 
l|Cutri. Rieke fc Lebofskvl Il984t ICondon et al Il99ll) . Mrk 
231 possesses a single compact nucleus surrounded by 
irregular 'rings' of recent star formation and a small tidal 
arm containing numerous blue starforming 'knots', and is 
therefore thought to be an advanced merger . The nucleus is 
unresolved in the radio l|Gondon et allll99ll) and at 11.7/im 
iMiles et all Il996l) . implying that the emission from the 
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Figure 1 - continued 



nucleus is powered by an AGN rather than a starburst. Mrk 
231 is classifie d as an AGN on the basis of both n ear-IR 
spectroscopy lllmanishi, Dudley fc Maloneyl 1200 if) and 
mid-IR spectroscopy (IRigopoulou et allll999l) . Conversely, 
the X-ray emission from this object, although weak, cannot 
be fitte d with a pure p ower law or thermal bremsstrahlung 
model (|lwasawa| ^99) . We compiled additional mid-IR 
data fro mlEigkg (EiZa) and additional su b- mm da ta from 
IRigopoulou. Lawrence fc Rowan-Robinsonl il99rJl and 



ICarico et al <198Sl) . We find that Mrk 231 contains both 
a luminous starburst and an AGN, where the starburst 
contributes ~ 70% of the 1 — 1000/mi IR luminosity. 
Interestingly, this result is in excellent agreement with in- 
dependent estimates of the starburst and A GN luminosities 
in M rk 231 derived from VLBI observations llLonsdale et all 
2003). From Figure [T] we can see that the AGN dominates 
the emission in the near- and mid-IR, in agreement with 
previous authors results, but that the starburst dominates 
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Figure 1 - continued 



at sub-mm wavelengths. 

13428+5608: This object is commonly known as Mrk 
273. It contains two nuclei with a single, long (41kpc) 
tidal tail and a 'ring' of star formation almost lOOkpc 
in diameter. The northern nucleus is resolved and is 
surrounded by unobscured star formation, whereas the 
south ern nucleus is redder and unresolved iScoville et all 
2000). The optical spectrum is that of a Sy2, and there 
are no broad lines in the near-IR spectrum, although 



there are strong narrow lines llVeilleux. Sanders fc Kiml 
1999). Radio continuum and HI 21cm observations of this 
source classify it as a pure starburst JCondon et allll99ll : 
ICarilli fc Tavlorl l2000l) . Converse ly , mid-IR spectral obser- 
vatio ns ilRlgOTaoTuouet a] 119991 : iLutz. Veilleux fc Genzell 
1999) classify it as an AGN. This object is a strong X-ray 
source, with a complex spectrum. The best fit is a 4 
component model, consisting of a direct and scattered 
AGN power law, a thermal component, and an iron line. 
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Figure 1 - continued 



Furthermore, the X-ray emission appears extended r ather 
than unresolved l|Levenson. Weaver fc HeckmanlEoOlh . We 
compiled an additional 350^im flux of 1004mJy ±230mJy, 
and an aditional 800/jm flux of 84mJy ±2 2m Jy, bo th from 
iRjgoDoulou. Lawrence fc Rowan-Robinsonl (|l99ft ). Our 
results are consistent with Mrk 273 containing both a star- 
burst and an AGN, in agreement with X-ray observations. 
The SED fit in Figure is consistent with a starburst 
interpretation in the near-IR, and also gives a good fit in 



the mid-IR although the fit predicts the mid-IR emission is 
due to a starburst rather than an AGN. 
14348-144^ : We compiled an additional 
800^m flux of 29mJy ± 10 mJy from 
iR.igoDoulou. Lawrence fc R.owan-Robinson| l)l99fift . This 
object is very luminous, where the luminosity is dominated 
by a starburst. Whilst this object does contain a luminous 
AGN, the AGN only contributes ~ 17% of the total infrared 
flux. 
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Figure 1 

15250+3609: This object has one bright central nucleus, 
and a second, much dimmer nucleus 0.7" from the centre. 
There is also a larg e ringlike structure 25kpc in diameter 
JScoville et all [20001. The optical spectrum is a composite 
of H II and LINER features llVeilleux. Sanders fc Kirr] 
Il999l) . iLutz. Veilleux fc Genzel lll999l) classify this object 
as a starburst based on its mid-IR spectrum. Our results 
however show that this object contains both a starburst 
and an AGN. This is one of only 3 objects in the sample 
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continued 

where the AGN provides more than half of the total IR 
emission. Interestingly, it can be seen from Figure Q that 
the sub-mm emission from this object is not dominated by 
a starburst, but instead arises in almost equal parts from 
a starburst and an AGN. We interpret this as arising from 
a heavily obscured AGN that is oriented nearly edge on to 
us, rather than an AGN with an extended torus. 
15327+2340: This object, commonly known as Arp220, 
is both the closest ULIRG to us and the most intensely 
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Figure 1 - continued 



studied. The optical spectru m was first thought to be 
Sy2-like ^Sanders et aj Tl988l. and this was interpreted 
as indicating an AGN power source. Conversely, the 
broad band optical-IR SEDs suggested a sta rburst as 
the main power source ij.Toseph fc WrighH I198ItI) . Obser- 
vations and modelling at other wavelengths revealed a 
more complex picture. Ground -based mid-IR spectroscopy 
l|Smith. Aitken fc Roch3ll989h suggested a hybrid nature 
for Arp220, including both starburst and AGN, whereas 



ISO mid-IR spectroscopy JSturm et all Il996h suggested a 
pure starburst. High resolution VLB A imaging dSmith et all 
1998) revealed a number of unresolved sources in the nuclei, 
inferred to be radio supernovae, supporting the presence of 
a starbu rst. Optical integral field observati ons of the central 
regions (lArribas. Colina fc Clementsll200lh showed that the 
spectrum is LINER-like in most regions, but w ith a Sy2-like 
spectrum in the brightest emission line region. iRiekel lll988h 
searched for hard X-rays in Arp220 with HEAO-A1 data, 
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Table 3. Ultraluminous Galaxies: Luminosities and Model Parameters 

Name x^ 3 L'j'g If* Lf£ N Age 6 SFR C ¥ MJ T? 0» 

L Q L Lq Myr A/Qyr -1 ° M kclvin 

00198- 7926 3~50 12.72^;^ 12.65+^3 < 11.91 64 - 71 130 ± 50 = 7.31+Jj;^ 36 ± 4 1.95 ±0.22 

00199- 7426 2.35 12.36^°;°^ 12.301°;°* 11.34+°;^ 26- 37 360 ± 50 8 - 48 to'.i2 27 ±7 1.96 ±0.40 
00262+4251 1.66 12.10l°'oi n - 64 iaoi n - 92 lo!o6 26 - 45 50 ± 30 < 15 8.42+. ,; 2 , , 25 ± 5 1.87 ±0.05 
00335-2732 1.33 12.01+°;°* 11.99+°;°| < 10.95 < 57 170 ± 20 7 - 07 -o'.02 42 ± 3 1-97 ±0.20 
01388-4618 2.08 12.11+.°;°| 12.031°;°^ < 11.35 > 26 200 ± 20 8.10±°-°;! 28 ± 4 1.92 ±0.10 
02364-4751 1.93 12.20+. , 01 12.15+. ;°,2 11.23+. ,;? 7 , > 10 260 ± 30 7 - 97 ±o.03 32 ± 7 1.91 ±0.20 
03068-5346 1.56 11.96±o.o2 n - 93 to.03 < 10.94 1.6-6.6 180 ± 40 7.62j:g;°J 35 ± 3 1.97 ±0.20 
04232+1436 1.41 12.07t°;^ 11.95+.°;:^ < 11.70 < 37 210 ± 30 7.87+JJ1J 31 ± 4 1.90 ±0.25 
05189-2524 2.20 12.16 + .jj;jj£ 11.99±q^| U.67±o H 6.6- 26 170 ± 30 25- 54 T.ei+jj; 2 ,® 33 ± 5 1.94 ±0.06 
06035-7102 1.00 12.24±°;°2 12.19+.°,;°,;!: n - 21 -ox>5 16- 26 290 ± 30 > 25 7.74^°;°^ 35 ± 2 1.95 ±0.03 
06206-6315 1.00 12.21+°;^ 12.18+°;"^ 11.05+.°;™ 16- 26 300 ± 30 < 45 7.93+.g;° 2 31 ± 2 1.92 ±0.02 
08572+3915 2.22 12.17^'^ 11.99j;^g 11.70+.°,;°J < 57 180 ± 20 7 - 04 to.06 45 ± 5 1-92 ±0.10 
09111-1007 1.92 12.0lto o2 ^-Oll ,; , 2 < 10.64 26- 37 192 ± 25 7.90to.ll 29 ± 4 1.94 ±0.28 
09320+6134 2.11 12.01+°;^ 11.98+.°;°i 10.86+.°,;°o > 16 160 ± 20 > 25 8.33+° °* 25 ± 3 1.93 ±0.03 
09583+4714 1.00 12.06 + .jj;jjg ll.gotg;" < 11.80 160 ± 70 34 ± 15 1.95 ±0.70 
10035+4852 1.00 11.98±2;^ U.9rtom < 11.50 180 ± 30 7 - 87 to'.i5 29 ± 5 1.95 ± 0.20 
10190+1322 1.00 12.10+° °2 12.01+°;°3 11.36l° °3 < 16 195 ± 15 7 - 85 -o.05 32 ± 3 1.97 ±0.15 
10494+4424 2.04 12.20+°^ 12.16+.°;°, 2 < H-06 < 26 280 ± 40 7 - 97 lo.05 32 ± 5 1.94 ±0.26 
10565+2448 1.23 Umt^l 12.03+.°;°,2 10- 8 Oto'.io 26-37 200 ± 20 7.85^2;^ 30 ± 2 1.95 ±0.20 
12112+0305 2.30 12.29+° °2 12.29+.°;° 2 < 9.50 16- 26 390 ± 20 7.94+°;°° 33 ± 2 1.95 ±0.10 
12540+5708 1.56 12.58; 1 ; ,; 04 12 - 4 2t°,:o6 l 2 - 05 to.oi 10- 16 470 ± 50 27- 54 7.95±°,;° 2 36 ± 2 1.94 ± 0.01 
13428+5608 1.99 12.15i2;°3 12.09+.°;°,3 11 - 21 -o'To 16- 26 260 ± 30 7 - 7 0±o'.02 31 ±4 1.95 ±0.08 
14348-1447 1.98 12.30+ ,; 01 12.18+. ,; 04 11.541°;^ 26 - 37 310 ± 30 < 32 7 - 97 t°,'.o2 32 ± 1 1.91 ±0.02 
14378-3651 1.00 12.12^0 02 12 -0 3 lo.oi < 11-36 16- 26 250 ± 40 7 - 81 to.20 34 ± 5 1.95 ±0.3 
15250+3609 2.32 12.08+° °4 11.681°;™ H- 8 6t°'°o 10- 37 70 ± 30 < 15 7.43+_°,;!j l | 37 ± 4 1.69 ±0.1 
15327+2340 2.01 12.14+°;°^ 12.14+. ,; , 2 < 9.90 26- 37 240 ± 30 8.32+° Q °l 27 ± 3 1.93 ± 0.02 
17132+5313 1.00 11.92+.°;°^ 11.841°;°^ 11.14+_° ) ;!j'J 10 - 26 130 ± 20 7 ^ 4 -o.05 31 ± 3 1.95 ± 0.03 
17208-0014 4.70 11.99inM 1L99 ±ao5 < 1100 > 16 150 ± 30 8 - 25 to'.25 26 ± 6 1.95 ±0.3 
18470+3233 1.00 la.OS+ij;^ 12.03+.°;°* < 11.42 6.6-37 200 ± 40 7.54+°!,!, 35 ± 6 1.95 ±0.05 
19254-7245 1.52 12.24t°-°o 12.17+.oo° < 11.60 170 ± 50 7 - 60 lo'.i4 34 ± 4 1.96 ±0.1 
19297-0406 1.00 12.4lt°;°g 12.39±°;°J < 11.50 < 57 480 ± 150 8 -l 5 ±o'.i5 31 ± 6 1.95 ±0.25 
19458+0944 1.73 12.40+°;°2 12. 25+°;°^ 11.83+°;°g ~ 26 280 ± 30 8 -"-o.i2 22 ± 2 1.89 ±0.01 
20046-0623 1.00 12.17+.o;jj4 12.11^°;^ < 11.60 < 26 210 ± 60 7 - 8 6t° °g 32 ± 9 1.96 ± 0.55 
20414-1651 1.18 12.22t° o2 i 2 -! 6 ! ,; ,! H-3ll° °o 6.6- 16 300 ± 15 7 - 49 ±o'.i4 38 ± 2 1.96 ±0.10 
20551-4250 1.95 12.05+°;° ) g 11.90+°;^ 11. 53+°;°^ 10- 57 160 ± 20 7 - 05 -o.05 44 + 3 1- 9 1 ±0.04 
21130-4446 1.65 12.14+.Q; , 2 12.11+.°°3 < 10.90 37-64 250 ± 15 7 -99lo.o3 30 ± 2 1.94 ±0.10 
21504-0628 1.00 12M±q%1 11.94+.°;°* 11.06t°;^ < 57 180 ± 20 7 -l 5 ±o'i5 42 ± 4 1.97 ±0.10 
22491-1808 2.38 12.16+° gi: 11.61+°;°^ 12.02+°,;°| ~ 26 60 ± 10 < 10 8.04+°;°^ 27 ± 3 1.80 ±0.02 
23128-5919 2.15 12.02+.°;° 2 11.93+.°;°,* ll- 30 to'.O7 26 - 37 150 ± 30 7 - 75 to.08 30 ± 3 1.93 ± 0.03 
23365+3604 2.90 12.19+2; 12.00+. ,; , 1 , ll- 73 to.o9 ~ 26 200 ± 20 < 10 7 - 8 ll"o.o3 33 ± 2 1.86 ±0.02 
23389-6139 1.22 12.16to.ol 12.15+. ,; , 1 , 10- 48 to.i8 16 - 26 290 ± 30 - 7 - 7 0±o'.09 33 ± 4 1.94 ±0.05 

Luminosities are the logarithm of the 1-1000/im luminosities obtained from the best fit combined starburst/AGN model, in units of 
bolometric solar luminosities. "Reduced x 2 of the combined SED fit b Starburst age c Star formation rate d Viewing angle of the AGN 
dust torus e Dust mass ■'Dust temperature 9 Emissivity index 



and concluded that an AGN similar to those in quasars 
could not be the power s ource . Recent observations with 
Chandra iClements et all 120021) detect weak hard X-ray 
emission from the nucleus, which has a possible origin in 
a weak AGN responsible for about 5% of the bolometric 
luminosity. An alternative explanation, with a much larger 
AGN contribution to the bolometric luminosity, is possible 
if the AGN is obscured by a Compton-thick screen of ~ 10 25 
cm -2 or more. There is some additional evidence for an 
obscured AGN i n Arp220, including sub-mm observations 
llHaas et"all200ll) and LWS spect ra that may indicate r > 1 
at ~ 100/im jFischer et al [l999l. We compiled add itional 
IR photometry for Arp220 from iKlaas et all l|l997ft . and 



additional sub-mm photometry from ICh ini et"al l)l98f^ : 
Rigopoulou. Lawrence fc Rowan- Robinson! j^9Bl and 
Durm^^^Ea1esn20MF ~Our results are consistent with a 
starburst interpretation for Arp220. The extensive IR data 
is well fitted by a pure starburst model, and the resulting 
upper limit on the IR luminosity of any AGN component 
is very strong. This upper limit on the AGN IR luminosity 
also allows us to set a lower limit in the UV optical depth 
of an IR luminous AGN in Arp220, of t uv = 1500. 
17208-0014'- This source has a single nucleus surrounded by 
a disturbed disk containing several compact star clusters, 
with a single tail. The optical spectru m is that of an HII 
region JVeilleux. Sanders fc Kimfll999ll and near-IR imag- 
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ing shows an exten ded nucleus, arguing against an AGN 
JScoville et allE oOO). This object lies in a crowded field, and 
high spatial resolution 100/im observatio ns with the Kuiper 
Airborne Observatory iZink et al 2000) measure a 100/im 
flux substantially below that from IRAS, implying that 
the IRAS fluxes for this obj ect are confused. The 100/im 
flux used in the fit is from IZink et a 3 (l2000M . This was 
the only object in our sample where the fit was relatively 
poor (Xred = 4.7), however it can be seen from Figure Q 
that this poorness of fit is due to inconsistencies between 
the model and the IRAS 60/im flux. As there are observed 
inconsistencies at 100/im we infer that this inconsistency 
at 60fim is also due to confusion. Overall, our results are 
consistent with a pure starburst interpretation for this 
object, although the AGN upper limit is weak. 
19254-7245: This object is more commonly known as the 
Super-Antena, and possesses a spectacular morphology. 
Two extremely long tail-like features extend on either side of 
the central regions, with a total length of 500kpc. These fea- 
tures can be explained as a result of a coplanar encounter 
betw een two massive disk galaxies l|Melnick fc Mirabell 
1990). There are two optical nuclei in the central regions, 
one with a Sy2 spectrum and the othe r with a starburst 
spectrum dMirabel . Lutz fc Mazal Il99ll). Further optical 
spectroscopy bv lCohm^LiparT^^ Vlacchettol <ll99ll) classify 
this source as an obscured AGN surrounded by starformhig 
clouds. The mid-I R spectrum of this obje ct classifies it 
as an A GN jLutz. Veilleux fc Genzell I1999T) X-ray obser- 
vations dPappa. Georeantopoulos fc Stewart |2000F) were 
inconclusive, but suggested an AGN was at least present 
in this object. We compiled for this object an extra flux a t 
6.7/im of 113mJy ±23mJy from lCharmandaris et all {2002). 
Of our sample the SED fit for this object is of reasonable 
quality, but deviates markedly from one point. We interpret 
this as being due to the very large angular size of this 
object, which is much larger than any other object in our 
sample. Whilst IRAS fluxes are likely to contain all the flux 
from the object, it is feasible that other observations have 
missed some flux. Our fit (Figure^ is a pure starburst, but 
the AGN upper limit is weak. 

224-91-1808: This system possesses a complex morphology, 
with at least 3 possible nuclei and large numbers of star- 
forming 'knots'. There are also two small bright tails, in 
which many of the st arforming knots are embedded. Based 
on the morphology, ICui et all ljl99lT) propose a multiple 
merger origin for this ULIRG. The mid-IR spectrum 
is consistent with a starbu rst l|Rigopoulou et al Il999t 
iLutz. Veilleux fc Genzell Il999l) . We however find that this 
object contains both a starburst and an AGN. This is one 
of only three objects where the AGN contributes more than 
half the total IR luminosity. In the mid-IR the starburst 
dominates, but with a significant AGN contribution. The 
contribution from the AGN in the sub-mm is comparable 
to the sub-mm emission from the starburst, implying that 
the AGN torus is being viewed almost edge on. 



6 DISCUSSION 

6.1 Limitations of the models 

Before discussing our results, we examine the suitability of 
the starburst and AGN models used in this paper for under- 
standing the nature of the power source in ULIRGs. Both 
sets of models are physical in nature, in that they do not 
assume an observational SED template for dusty star for- 
mation or AGN activity, and instead produce a predicted 
spectrum based on radiative transfer calculations of opti- 
cal/UV light from a starburst or AGN propagating through 
a dusty medium. They are therefore well suited to investi- 
gating the IR emission from a population whose properties 
may vary markedly between individual objects. A complete 
d escription of the starburst and AGN models can be f ound 
hrjEfstatluo^ ^ fcSiebenmorgenl i2000fl and 

lEfstathiou fc Rowan-Robinsonl jl995l) respectively. 

It is however important to note that these models have 
drawbacks. In particular, the models cannot account for sub- 
mm emission from an extended (hundreds of parsecs in- 
stead of tens of parsecs) dust torus surrounding an AGN. 
If such a torus did exist in our sample the sub-mm emission 
from the torus could be mistakenly ascribed to star forma- 
tion, leading to a substantially overestimated star formation 
rate. There do exist radiative transfer models for cool, ex- 
tended dust surrounding an AGN. To test whether all or 
part of the sub-mm emission from our sample could arise 
from such a toru s , we fitted the extended torus models of 
iRowan-Robinsonl jl995l) to the objects in our sample, both 
on their own and in combination with the starburst mod- 
els described previously. Although in all cases the extended 
torus models were clearly rejected, we cannot completely dis- 
count this possibility, as the extended torus models do not 
cover a wide range of torus geometries and lines of sight. 

Despite their limitations, the models appear to work 
very well. In a sample of 41 objects the models produce 
a good fit in all but one case. As described in Section |S] in 
this case there is reason to believe that the relative poorness 
of fit may be partially due to the quality of the data. The 
models successfully reproduce both a wide range of ULIRG 
SED shapes over the wavelength range 0.4 — 1000/im (if the 
optical data are used as upper limits), and also reproduce 
spectral features such as the 7.7/im PAH line, and the 9.7/im 
silicate absorption feature. We therefore find that, under the 
caveats described previously, these models are well suited to 
studying the dust-shrouded power source in ULIRGs. 

6.2 Starburst and AGN activity in ULIRGs 

6.2.1 What powers ultraluminous infrared galaxies? 

Although ULIRGs have been studied for nearly two decades, 
the nature of the power source behind the IR emission has 
remained controversial. The most recent results on the power 
source in ULIRGs have come from mid-IR spectroscopy, 
which show tha t most ULIRGs [~ 80% ) are powered mainl y 
by starbursts iRigopoulou et al Il999t iGenzel et all 1199s!) . 
but that at least half of local ULIRGs show evidence for 
both starburst and AGN activity. The AGN fractional lu- 
minosity as a function of total IR luminosity has been 
studied si nce IRAS galaxies w ere discovered. Optical spec- 
troscopy JVeilleux et all Il995l) found that the fraction of 
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Figure 2. AGN luminosity in the IR plotted against starburst luminosity in the IR. The ULIRGs from our sample are plotted as 
points, with a downward arrow indicating an upper limit on the AGN luminosity. Triang les are PG Q SOs, and crosses are those galaxies 
with detections in all four I RAS bands, both the se samples are taken from iRowan- Robinson! 120001) . Diamonds are the Hyperluminous 
Infrared Galaxies taken from lFarrah et a 3 1200211 . 



IRAS sources with AGN spectra, and the fraction of Seyfert 
galaxies amongst the AGN increases with increasing IR lu- 
minosity, reaching values of 62% and 54% respectively at IR 
luminosities > 10 12 Lp>. Near-I R spectroscopy of ULIRGs 
JVeilleux. Sanders fc Kimlll999l) shows that the fraction of 
ULIRGs with signs of AGN activity is at least 20% - 25% 
but rises abruptly to 35% — 50% for objects with Li r > 
10 12 ' 3 Lp,. Recent ISO sp ectroscopy of a small sample of 
ULIRGs jTran et all 1200 ll) found half of the sample to be 
starburst dominated and half to be AGN dominated. They 
also showed that, at IR luminosities below lO 12 ' 4 !/©, most 
ULIRGs were starburst dominated, with the starburst con- 
tributing around 85% to the IR emission. At IR luminosi- 
ties above W 12 4 Lq the AGN contribution was much higher, 
with the starburst contributing about 50% of the IR emis- 



sion on average. Starburst dominated systems were found 
up to luminosities of around W 12 65 Lq. 

From our results, we find that both starburst and AGN 
activity are central in understanding the properties of local 
ULIRGs. All 41 objects in our sample contain a very lumi- 
nous starburst that contributes significantly (at least 28%) 
to the total IR luminosity. There are no purely AGN pow- 
ered systems in the sample. 23/41 objects have measured 
AGN luminosities where the AGN contributes significantly 
to the total IR emission, the remaining 18 objects have only 
upper limits on the AGN contribution. In three cases the 
AGN supplies more than half of the total IR emission, in 
the remaining cases the starburst is the dominant contrib- 
utor. We therefore find that previous estimates of the frac- 
tion of local ULIRGs that are starburst dominated based on 
mid-IR spectroscopy have underestimated the true fraction, 
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Figure 3. Total IR luminosity plotted against (left) starburst l uminosity and (right) AGN luminosit y. Our ULIRGs are plotted as 
points, triangles are PG QSOs taken from lRowan- R.obinsonl l200Ct> . and diamonds are the HLIRGs from lFarrah et all l200 2bl 



and that this fraction is ~ 90% rather than 80%. The mean 
starburst fractional luminosity for the sample is 82%, span- 
ning the range 28% to ~ 100%. Overall, the IR emission 
from ULIRGs as a class is either starburst in origin with a 
negligible AGN contribution, or arises from a combination 
of starburst and AGN activity with the starburst usually 
contributing the largest fraction. Given our sample size, the 
fraction of purely AGN-powered local ULIRGs must be less 
than ~ 2%. The derived star formation rates range from 
5OM0yr _1 to 5OOM0yr~ and the derived dust masses span 
the range 10 7 < Mq < 10 9 . These ranges are at least 1 order 
of magnitude higher than those observed in normal galaxies, 
indicating that all ULIRGs are going through a major star 
forming episode. It proved difficult to constrain the line of 
sight to the AGN torus, with meaningful constraints only 
achieved for a few objects. Although these values are dis- 
cussed in section |K| we cannot draw any general conclusions 
about relative AGN orientation in ULIRGs. 

Figure [5] shows a plot of starburst luminosity vs. AGN 
luminosity. Also plotted are the starburst and AGN lumi- 
nosities of 10 Hyperlu minous Infra r ed Gal axies (HLIRGs, 
diamonds) taken from iFarrah et all l|2002b|) . as well as 22 
PG QSOs (triangles) and 78 galaxies with detections in all 
four IRAS bands (crosses). Th e two latter samples are taken 
from iRowan-Robinsonl 820001) and rescaled to Ho = 65 km 
s^ 1 Mpc . The figure shows that starburst and AGN lu- 
minosities are correlated, albeit with a large scatter, over 6 
orders of magnitude in IR luminosity and over a wide range 
of galaxy types. Therefore, there may be a common physi- 
cal factor governing the IR luminosity of both starbursts and 
AGN in these galaxies. We postulate that this factor is the 
available quantities of gas and dust in the nuclear regions of 
these systems, as this affects both the rate and duration of 
black hole accretion, and the number of stars that can form. 
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Figure 4. Total luminosity in the IR plotted against AGN frac- 
tion. Key to the symbols is the same as for Figure 131 



We next examine trends in starburst and AGN lumi- 
nosity against total IR luminosity. Figure |3] shows total IR 
luminosity plotted against starburst luminosity and AGN 
luminosity. Also plotted are the PG QSO and HLIRG sam- 
ples, however in this section we only discuss the ULIRGs 
plotted in these figures; comparisons between the ULIRGs, 
PG QSOs and HLIRGs can be found in ffljX2l The left hand 
panel shows that the ULIRG starburst luminosity is strongly 
correlated with the total luminosity. From the right-hand 
plot in Figure[3]we can see that, although there is a trend for 
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more luminous ULIRGs to contain a more luminous AGN, 
the trend is less clear than for the starburst luminosity. Fig- 
ure 0] shows total IR luminosity plotted against fractional 
AGN luminosity. It can be seen that there is no trend for 
increasing AGN dominance with increasing total luminosity 
in ULIRGs, contrary to previous claims. We find that these 
claims were based on finding generally more luminous AGN 
in more luminous ULIRGs, but that there is no evolution 
of fractional AGN luminosity with total luminosity in local 
ULIRGs. 

6.2.2 Multiple starbursts and multiple mergers 

The number of individual starburst events that occur over 
the lifetime of a ULIRG is an important parameter in under- 
standing starburst triggering and evolution in all classes of 
active galaxy. Similarly, the number of progenitors in a typ- 
ical ULIRG merger can be used to estimate which fraction 
of ULIRGs occur in the field galaxy population between 2 or 
3 progenitors, and which fraction occur in compact groups 
of galaxies between multiple progenitors. 

Both these parameters have proved difficult to estimate. 
Numerical simulations of mergers between two galaxies have 
shown that multiple, discrete starburst events can be trig- 
gered even in double mergers. The timing of the starbursts is 
thought to depend on the progenitor galaxy morphologies, 
with early starbursts (before the galaxies start to merge) 
occurring in bulgeless disk galaxies, and late starbursts (af- 
ter the galaxies have started to merge) occurring in disk 
galax ies with a bulge component iMihos fc Hernauisdll994l 
1996). Simulations of multiple galaxy mergers have how- 
ever shown that more than two starburst or AGN nuclei 
is evi dence that there were mor e than two merger progen- 
itors (iTaniguchi fc Shioval Il99sh . and that re petitive star - 
bursts are characteristic of multiple mergers feekkil 12001). 
Obs ervationally, multiple mergers have been linked to Arp 
220 JPiamond et a]||l989Tl and to lar ger samples of ULIRGs 
jBorne et all2000l:lFarrah et al200lh . Other authors have es- 
timated the lif e time o f a single starburst event in ULIRGs. 
iThornlev et all <|2000T) derive a starburst lifetime range of 
10 6 — 10 7 years, and infer that galactic superwinds pro- 
duce d by supernovae ma y be responsible for the short dura- 
tion. iGenzel et all il99ST> derive an upper limit on the star- 
bust lifetime in ULIRGs of ~ 10 s years. The lifetime of 
a ULIRG has been previously estimated to lie in the range 
10 s -10 9 years jMurphv et al200ll:lfVrah et al200lh based 
on the observed range of morphologies in large samples of 
ULIRGs, the apparent discrepancy between the ULIRG and 
starburst lifetimes has led to the suggestion that multiple 
starbursts and mul tiple mergers may be common in ULIRGs 
jFarrah et a1l200ll) . 

The starburst ages derived for our sample are listed in 
Table[3] As our sample is large and the selection is robust we 
are confident that this range in ages is representative of the 
ULIRG population as a whole. Excluding those objects with 
upper or lower limits on the starburst age, the derived ages 
span the range 1.6 x 10 6 ^ years sC 7.1 x 10 7 , with most lying 
in the range 1.0 x 10 7 < years < 3.7 x 10 7 . These lifetime 
ranges are clearly inconsistent with a single starburst event 
powering the IR emission for the lifetime of a ULIRG. From 
section IB . 2 . 1 1 however . starbursts are present in at least 98% 
of ULIRGs and in most cases dominate the IR emission; 



AGN activity cannot therefore explain this discrepancy in 
lifetimes. We therefore find that multiple starburst events 
must be common in the ULIRG population, and by inference 
that multiple mergers must be common as well. 

From this analysis, coeval starbursts with different 
ages are feasible in ULIRGs, and it is important to note 
the effects this may have on our SED fitting approach. 
Whilst the available photometry is sufficient to separate 
starburst and AGN emission, it is not sufficient to resolve 
the starburst emission into multiple components. The de- 
rived age ranges in Table should however encompass the 
age ranges of any coeval starbursts, due to the effect of 
age on the shape of the starburst SED. As the age of 
the starburst increases then (with reference to figur e 3 of 
lEfstathiou. Rowan-Robinson fc Siebenmorgenl l|2QQQ,)), the 
peak of the SED shifts to longer wavelengths, PAH features 
become stronger, the 9.8pm silicate absorption feature be- 
comes shallower, and the optical/UV flux increases. Hence, 
the mid- and far-IR photometry, in constraining the strength 
of the mid-IR spectral features and the peak of the SED re- 
spectively, provide a lower limit to ages of any coeval star- 
bursts, whereas the optical photometry provides an upper 
limit. Therefore, even if some or all of the ULIRGs in our 
sample contained two or more starbursts with different ages, 
we are confident that the starburst age ranges given in Table 
[3] encompass the ages of these starbursts. 

6.3 Observational diagnostics of ULIRGs 

6.3.1 Mid-infrared spectroscopy 

Several groups have constructed diagnostics for the power 
source in ULIRGs based on mid-IR spectra from ISO ob- 
servations, specifically using the relative strengths of the 
Unitentified Infrared Emission Band (UIB) features lying 
between 5(j,m and 12^im. Currently, t he most prev a lent o f 
these diagnostics is that proposed by IGenzel et a 3 Jl998l) . 
In this diagnostic the ratio of the 7.7/im emission, thought 
to arise from PAHs, relative to the continuum emission at 
the same wavelength is used to differentiate between star- 
burst and AGN emission. A high value of F7.7/C7.7 (> 1) 
indicates starburst power whereas a lower value of -F7.7/C7.7 
indicates AGN power. This is generally borne out by mid-IR 
obse rvations of lower luminos ity starburst and AGN systems 
fe.g. iRjgopoulou et all (|l999lh 'l. There are however potential 
problems with this approach. Firstly, since the active regions 
in ULIRGs are small and dusty it is not clear that even 
mid-IR observations can penetrate to the central regions, 
although it is argued that this is unlikely to be significan t 
iGenzel et allll998l: IRjgopoulou et all999t iTran et a1l200lTl . 
Secondly, some metal-poor dwarf starburst galaxies show no 
UIB features, but as ULIRGs are not metal poor thi s is also 
argued not to be important (Rigopoul ou et al ll999). 

By comparing the results from diagnostics using the 
F7.7/C7.7 ratio to o ur results for the 22 ULI RGs common to 
our sample and to Rigo poulou et all il999t> . we can exam- 
ine the reliability of the F7.7/C7.7 technique. In Figure|5]we 
plot F7.7/C7.7 against starburst luminosity and fractional 
AGN luminosity. From the left hand plot it can be seen 
there is, at best, a weak trend for more luminous starbursts 
to have lower values of F7.7/C7.7. The immediate possibil- 
ity is that the mid-IR diagnostics are correct and our SED 
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fitting approach is wrong, however our SED fits reproduce 
the observed F7.7/C7.7 ratios in all 22 objects, as well as the 
observed SED from the near-UV to the sub-mm. We there- 
fore infer that, of the two approaches, ours is more robust. 
Therefore, it seems that the relative strength of the 7.7/im 
PAH feature is not a reliable indicator of the luminosity 
of the starburst, but merely of its presence. We postulate 
that the discrepancy between the nature of the 7.7/im PAH 
feature in ULIRGs and in less luminous starburst galaxies 
is due to two reasons. Firstly, in ULIRGs it is more likely 
than in moderate starbursts that the obscuration is so high 
that spectral features are not apparent in the mid-IR. Sec- 
ondly, the intense UV radiation in ULIRG-like starbursts 
may be capable of dissociating the hydrocarbon bonds in 
PAHs. Both these reasons would explain why more lumi- 
nous starbursts have smaller F7.7/C7.7 ratios. 

More serious however is that the right hand panel of 
Figure |S] shows no correlation between -F7.7/C7.7 and the 
fractional AGN luminosity, from which we infer that the 
F7.7/C7.7 ratio on its own is not an accurate indicator of 
the overall power source in U LIRGs. This discrepancy be- 
tween our results and those of IRigopoulou et all il99Sl) can 
also be explained by the more extreme nature of ULIRGs. 
From our results, all ULIRGs contain a starburst, explain- 
ing the prevalence of PAH features in their mid-IR spec- 
tra, however the majority of ULIRGs are composite systems 
and may contain a heavily obscured starburst or AGN that 
does not emit significantly in the mid-IR but that does emit 
strongly in the far-IR/sub-mm. As an example we consider 
IRAS 00262+4251. This object, with F7.7/C7.7 = 4 063 i s 
classified as a pure starburst bv IRigopoulou et all il998t) . 
From our SED fit however (Figure EJ, the AGN does not 
contribute in the mid-IR as the torus is viewed almost edge 
on and the corresponding obscuration in the mid-IR is very 



high. At wavelengths between 50/mi and ~ 120^im the AGN 
contribution is however dominant. Although this object does 
contain a starburst, overall the emission is mainly due to 
AGN activity. We conclude that the F7.7/C7.7 ratio on its 
own is only an indicator of whether a ULIRG contains a 
moderately obscured and moderately luminous starburst. It 
cannot be used to probe either heavily obscured or very 
luminous starburst activity, or AGN that are either isotrop- 
ically obscured or oriented nearly edge on relative to us. As 
such, it cannot be used on its own to probe the overall power 
source in ULIRGs. 



6.3.2 The radio-infrared correlation 

There is observed to exist a tight correlation between 
radio and infrared flux for both normal and active 
galaxies, extending over several ord ers of magnitude 
llHelou. Soifer fe Rowan-Robinsonlll985l) . The physical ori- 
gin of this relation is believed to be a population of rel- 
ativistic electrons accelerated by supernov a remnanats pro- 
duced in regions of massive star formation (lHarwit fc Pacini 
Il975l) . This relation is usually expressed using the param- 
eter q, the logarithm of the ratio of infrared flux to radio 
flux. For starburst galaxies and ULIRGs, q is obse rved to 
lie in the range 2.0 < q < 2.6 iCondon et alll99lT> . Radio 
Loud Quasars on the ot her hand have ty pical q values in 
the range < q < 1 (e.g. iRov et al dl998l) ). whereas Radio 
Quiet Quasars have q values lying in a similar range to loca l 
starburst galaxies and ULIRGs dSopp fc Alexanderlll99ll) . 
It has therefore been suggested that the radio-IR correla- 
tion in both U LIRGs and Radio Quiet Quasars is due to 
star formation iCondon et alll99ll : ICram. North fc Savage! 
Il992l) . 

Using the results from this paper we can examine the 
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(Left) Total IR luminosity plotted against total radio luminosity. (Right) Total IR luminosity vs. q parameter. Open squares 
crosses are q 3 tarburst> an d small filled squares are qAGN- A downward arrow indicates an upper limit. 



Table 4. Ultraluminous Infrared Galaxies detected in the radio 



Name 


FlAGHz 


qtotal 


Qstarbur st 


qAGN 


00262+4251 


28.0 


1.92 


1.32 


1.80 


00335-2732 


11.25 


2.60 


2.58 


< 0.70 


04232+1436 


29.39 


2.10 


2.10 


< -0.4 


05189-2524 


28.8 


2.66 


2.62 


1.53 


08572+3915 


4.3 


3.21 


3.16 


2.27 


09320+6134 


170.1 


1.85 


1.85 


0.00 


09583+4714 


36.5 


1.90 


1.87 


< 0.60 


10035+4852 


28.1 


2.26 


2.26 


< -0.40 


10494+4424 


21.2 


2.28 


2.28 


< -0.30 


10190+1322 


16.8 


2.38 


2.37 


0.48 


10565+2448 


56.9 


2.38 


2.37 


0.45 


12112+0305 


23.3 


2.57 


2.57 


< -0.37 


12540+5708 


308 


2.03 


2.00 


0.88 


13428+5608 


144 


2.22 


2.22 


-0.15 


14348-1447 


35.8 


2.26 


2.17 


1.53 


14378-3651 


33.8 


2.29 


2.29 


< -0.50 


15250+3609 


14.5 


2.70 


2.21 


2.53 


15327+2340 


326.5 


2.42 


2.42 


< -1.0 


17132+5313 


29.3 


2.34 


2.32 


0.98 


17208-0014 


81.8 


2.07 


2.07 


< 0.76 


18470+3233 


12.4 


2.50 


2.50 


< -0.08 


19458+0944 


14.0 


2.52 


2.37 


1.99 


20046-0623 


14.6 


2.42 


2.42 


< -0.16 


20414-1651 


23.4 


2.32 


2.32 


-0.23 


21504-0628 


13.8 


2.41 


2.40 


0.62 


22491-1808 


5.9 


2.87 


2.33 


2.72 


23365+3604 


27.2 


2.47 


2.31 


1.93 



1.4GHz radio fluxes were taken from the NVSS catalogues and 
are in mjy. The q values were calculated using Equation 151 

origin of the radio-IR correlation for ULIRGs, and the fac- 
tors that may bias it. 1.4GHz fluxes were obtained for most 
of the objects in our sample from the NVSS catalogues, these 
fluxes are presented in Table|l| The left hand panel of Figure 



HJ shows total IR luminosity plotted against 1.4GHz lumi- 
nosity, and a weak, though positive, correlation can be seen. 
Figure [7] shows the starburst and AGN luminosities plotted 
against 1.4GHz luminosity. The correlation between the ra- 
dio luminosity and the starburst luminosity is stronger than 
for the total luminosity, whereas no correlation is observed 
with the AGN luminosity. We thus confirm that the radio-IR 
correlation in ULIRGs is due to star formation. 

We have further investigated this relation using the q 
parameter. To allow future comparisons with galaxy samples 
which do not benefit from accurate estimates of IR luminosi- 
ties from SED fits, we express q as: 

.-««(£) m 

where F1.4 is the 1.4GHz flux and Feo is the rest-frame 60/im 
flux. Rest-frame 60/im fluxes were extracted from the best- 
fit total, starburst and AGN SEDs. In Table 0] we present 
values of q calculated using these fluxes. We derive (qtotal) = 
2.36 ± 0.06, (qstarburst) = 2.29 ± 0.06 and (q A G N ) = 1.29 ± 
0.25, which are in agreement with previous estimates for 
ULIRGs. Individual q values are plotted against total IR 
luminosity in the right hand panel of Figure |5] It is evident 
that qtotal and q s tarburst have only a small scatter about a 
mean value, but that the values of qAGN appear randomly 
distributed. 

If the radio-IR correlation in ULIRGs is due to star 
formation, then what is the origin of the scatter about the 
mean values of q s tarburst and qtotal'? Broadly, there are four 
possibilities. Firstly, the IMF of the starburst may be skewed 
towards or away from producing high-mass stars, thus pro- 
ducing an over- or underabundance of radio supernova rem- 
nants, thus causing a scatter. This possibility is impossible 
to investigate within the context of this paper. Secondly, 
there may be contamination of the 1.4GHz luminosity from 
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Figure 7. (Left) Starburst luminosity plotted against 1.4GHz radio luminosity. (Right) AGN luminosity plotted against 1.4GHz radio 
luminosity. 



an obscured AGN, leading to an artificially suppressed q 
value. We argue that this possibility, although undoubtedly 
present in the ULIRG population, is not the main cause of 
the observed scatter. If there were a significant number of 
obscured AGN causing the scatter, then we would expect 
to see some correlation between AGN luminosity and radio 
luminosity, and yet none is observed. Even if only the most 
IR- luminous AGN are plotted in Figure |S| no such correla- 
tion appears. The third possibility is that the starburst may 
be too young to have formed a significant population of radio 
supernova remnants, thus causing an artificial amplification 
of the q parameter. We argue that this is a negligible ef- 
fect in the ULIRG population for two reasons. Firstly, most 
of the starburst ages lie in the range lOMyr - 37Myr, by 
which time a significant population of supernova remnants 
will certainly have formed. Secondly, if the youngest star- 
bursts are removed from the left-hand plot of Figure [7| then 
the scatter becomes bigger rather than smaller. The fourth 
possibility is an old relativistic electron population from a 
previous starburst event. As described in section 16.2.21 it 
is likely that ULIRGs undergo multiple, discrete starburst 
events with lifetimes in the range 10 6 — 10 8 years. Although 
radio supernova remnants have lifetimes of only 10 3 years at 
best, the relativistic electrons they produce have lifetimes of 
the order 10 years (ICondonl H992). This lifetime lies at the 
upper end of the lifetime range of a single starburst event, 
and is comparable to the total lifetime of a ULIRG. We thus 
conclude that most of the scatter around the radio-IR cor- 
relation for ULIRGs is due to skewed IMFs in the starburst, 
and relativistic electrons from a previous, separate starburst 
event. 



6.4 ULIRG evolution 

6.4.1 Do ULIRGs evolve into QSOs? 

The fir st picture of ULIRG evolution to be proposed was 
that of ISanders et al il988fl . which asserted that ULIRGs 
as a class evolve into optical QSOs. According to this pic- 
ture (hereafter referred to as the S88 picture), interactions 
and mergers between gas rich spirals transport gas to the 
central regions of the galaxies. This central gas concentra- 
tion fuels starburst activity, and in the latter stages of the 
merger commences the fuelling of a central supermassive 
black hole, which rapidly comes to dominate the IR luminos- 
ity of the ULIRG. In the last stages the dust shrouding the 
black hole is blown away and the ULIRG evolves into an op- 
tic al QSO. A qualitative l y simi lar picture was also proposed 
bv lCanalizo fc Stockton! I2001I). A more recent ev olutionary 
scenario for ULIRGs is that of Far rah et all (1200 ll) . also later 
suggested bv lTacconi et aj i2002l) . In this scenario ULIRGs 
are not simply the dust shrouded precursors to optical QSOs 
but instead are a diverse population with a broad range of 
properties and evolutionary paths. Although some ULIRGs 
would still evolve into QSOs, this number would be a small 
and non-representative subset of the ULIRG population. 

Previous studies of the evolution of ULIRGs have high- 
lighted some problems for the S88 picture. In the mid-IR , 
spectroscopic stud ies of ULIRGs iRigopoulou et all Il999l : 
iGenzel et all Il998l) find no evidence that the advanced 
merger systems are more AGN-like, based on the ULIRG nu- 
clear separations. That AGN activity becomes more preva- 
lent with merger advancement is a natural prediction of 
the S88 picture. Furthermore, a recent spectro scopic study 
of th e stellar and gas kinematic in ULIRGs jGenzel et all 
l200ll) found that ULIRGs will likely evolve into only ~ 
L* ellipticals, and not into the most massive ellipticals 
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Figure 8. ^25/^60 colour plotted against (left) AGN luminosity, and (right) fractional AGN luminosity. 



seen locally. A corollary to this are the results from a 
large scale imaging and spectr oscopy survey of ULIRGs 
JVeilleux. Kim, fc Sandersl2002l) . in which it was found that 
nearly half the sample may not evolve into an optical, post 
ULIRG AGN, and that merger induced QSO activity might 
only take place if both merger progenitors had L* luminosi- 
ties or greater. 

Using the results from this paper, and from previous 
authors, we can examine these two evolutionary scenar- 
ios. We first consider the growth of Supermassive Black 
Holes (SMBHs) in ULIRG mergers, and relate this to the 
SMBH masses seen in QSOs, and to the level of AGN ac- 
tivity seen in our sample. Nearly all optical QSOs with 
measured SMBH masses have Mbh ^ 10 s Mp, in some 
cases reaching Mbh ~ 10 10 M© jMcLure fe DunlodEoOlt 
iGu. Cao fc Jiandl200ih . hence it seems likely that a SMBH 
in this mass range is necessary for QSO activity. The 
growth of SMBHs in ULIRG merg ers has been studied by 
iTaniguchi. Ikeuchi fc Shioval dl99St) . They conclude that a 
SMBH of mass ~ 10 8 Mq can form if any of the progenitor 
galaxies contains a 'seed' SMBH of mass > 1O 7 M under- 
going efficient Bondi type gas accretion over ~ 10 8 years. 

Hence, if the S88 picture is correct and ULIRGs as a 
class evolve into QSOs, then the range of SMBH masses seen 
in spiral galaxies (10 6 < M B h(M q ) < 10 8 ) would have to be 
transformed by the merger into the SMBH mass range seen 
in QSOs (10 8 < Mbh(M q ) < 10 10 ). This would require 
very efficient accretion onto the black hole for the entire 
duration of the merger. This itself implies that the lifetime 
of the IR-luminous AGN phase will be significantly greater 
than the lifetime of any single starburst event. Furthermore, 
in order to achieve the SMBH mass range seen in QSOs, 
then the AGN in a ULIRG observed at a random point 
during the merger would in many cases harbour a SMBH of 
mass Mbh Z 1O 8 M0. Overall therefore, if the S88 picture 



is correct, then most ULIRGs should contain a long-lived 
AGN which, by virtue of efficient accretion onto a SMBH 
with mass comparable to those seen in QSOs, should have 
an IR luminosity comparable to the bolometric luminosity 
seen in QSOs. Our results however do not support this, as 
we find all ULIRGs to contain a very luminous starburst, 
with only about half containing a luminous AGN, and that 
less than 10% of ULIRGs are AGN dominated. Whilst we 
cannot rule out the S88 picture based on AGN activity alone, 
it seems unlikely that the level of AGN activity seen in our 
sample could transform the BH mass range in spirals into 
that seen in QSOs. 

We next consider the nature of the so-called 'warm' 
ULIRGs, those with F 25 /F 60 > 0.25, (where F 25 and F 60 
are the IRAS 25/rni and 60/im fluxes respectively). Objects 
with 'warm' infrared colours have been previously shown 
to be more likely to contain IR-luminous AGN. Hence, if 
the S88 picture is true, then the greater prevalence of AGN 
in 'warm' ULIRGs than in cool ones means that 'warm' 
ULIRGs will be, on average, nearer the end of the ULIRG 
phase than 'cool' ULIRGs, and closer to evolving into opti- 
cal QSOs. Therefore, 'warm' ULIRGs should have a higher 
fractional AGN luminosity on average than 'cool' ULIRGs 
as the AGN comes to dominate the bolometric emission and 
the starburst dies away. Figure |H| shows Fzs/Fso colour plot- 
ted against AGN luminosity and fractional AGN luminosity 
for the objects in our sample. Values of Fzs/Fqo are pre- 
sented in Table As expected, those objects with a higher 
value of F25/F60 are more likely to contain a more lumi- 
nous AGN. There is however no trend of increasing AGN 
fractional luminosity with increasing IR colour. 

Overall therefore, our results contradict the S88 picture, 
and imply that ULIRGs as a class do not evolve to become 
QSOs. Instead, the wide range of starburst and AGN lu- 
minosities derived for our sample imply multiple possible 
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patterns of starburst and AGN activity over t he lifetime of 
the ULIRG, in line with the lFarrah et all (1200 ll) picture. The 
overall level of AGN activity is consistent with sufficient ac- 
cretion to transform the range of black hole masses seen in 
spirals into the range of black hole masses seen in ellipticals, 
without recourse to an optical QSO phase. As a small num- 
ber of our sample are AGN dominated it is likely that some 
ULIRGs do become optical QSOs, but the small number of 
such systems in our sample argues that such ULIRGs are 
rare, and thus unrepresentative of the ULIRG population as 
a whole. 

6.4-2 ULIRG evolution with redshift 

ULIRGs, despite being rare in the local Universe, are a 
cosmologically important galaxy population. At low red- 
shifts the ULIRG luminosity function shows stro ng evolu- 
tion with redshift dVeilleux. Sanders fc Kimlll999l) . In addi- 
tion, deep sub-mm surveys find that systems with ULIRG- 
like luminosities domi nate the global ene r gy density of the 
Universe at z ^ 1 jBareer et aJ Il998l : iLillv et all Il999t 
IScott et ail2002t iFox et alH2002T) . Determining how systems 
with ULIRG-like luminosities evolve with redshift is there- 
fore necessary to measure, for example, the overall star 
formation histo ry of the Universe , currently a hotly de- 
bated topic (e.g. (Madau et al lll99d) :lR owan-Robinson et all 
J1997M . 

It is however not known whether the IR luminosity of 
these high redshift sources arises due to mergers between 
two or more large spiral galaxies, as in local ULIRGs, or via 
a different mechanism. If the former is true then the strong 
evolution seen in the ULIRG luminosity function will most 
likely be due to mergers between two or more large spiral 
galaxies becoming more common with increasing redshift. 
If the latter is true then this indicates that the strong evo- 
lution seen in the ULIRG luminosity function is probably 
due to different galaxy formation processes becoming more 
important at high redshifts. Currently, the only sample of 
IR luminous sources at high redshift with resolved starburst 
an d AGN components is the sample of HLIRGs presented 
bv lFarrah eTaf J2002rj) . On morphological grounds there is 
some evidence that the trigger for the IR emission in local 
ULIRGs and high-z HLIRGs is the same; nearly all UL IRGs 
show signs of mor phological disturbance feorne et alll200ol 
Farr ah et "a1|200lD as do approximately half of HLIRGs 
jFarrah et al l2002a s |. 

In section 16.2.11 we showed that the broad correlation 
between starburst and AGN luminosities for a wide range 
of galaxy types, including ULIRGs and HLIRGs, was most 
plausibly explained by the starburst and AGN luminosi- 
ties being governed by the gas and dust masses in the nu- 
clear regions. We now examine whether the evolution of 
the starburst and AGN components are the same or dif- 
ferent between ULIRGs and HLIRGs, to determine whether 
HLIRGs are the high redshift analogues of ULIRGs or a dif- 
ferent galaxy population, by co mparing the ULIRG s from 
our study with the HLIRGs from lFarrah et all (l2002bl) . Both 
samples span a similar range in fractional AGN luminos- 
ity, ranging from almost pure starbursts to AGN dominated 
systems. Amongst the ULIRGs most systems are starburst 
dominated with a mean starburst fraction of 82%, whereas 
in the HLIRG population approximately half the systems are 



AGN dominated, with a mean starburst fractional luminos- 
ity of ~ 35%. This difference could be because a more lumi- 
nous AGN is required to generate the higher total luminosi- 
ties in HLIRGs, however starburst dominated HLIRGs are 
found amongst the most luminous members of the HLIRG 
population. This then constitutes a distinct difference be- 
tween the two samples. Amongst local ULIRGs, starbursts 
are the dominant contributor to the total luminosity whereas 
amongst HLIRGs, although a starburst could in principle 
accomplish the same thing, this is not observed and AGN 
activity is much more prevalent than in ULIRGs. This can 
also be seen in Figures |3] and 0] in which the starburst and 
AGN luminosities are compared to the total IR luminosities 
in ULIRGs, HLIRGs and PG QSOs. Furthermore, it can also 
be seen from these figures, particularly the right-hand panel 
of Figure |S| that the pattern of starburst and AGN activity 
in HLIRGs is much more reminiscent of the pattern seen in 
QSOs than in ULIRGs. 

Based on this comparison, we conclude that IR- 
luminous galaxies at z ~ and z ^ 1 are physically dif- 
ferent galaxy populations. We speculate that this difference 
is due to different galaxy formation processes in the low and 
high redshift Universe, and that the trigger for starburst 
and AGN activity in high redshift IR-luminous sources is 
more similar to the trigger for QSO activity, rather than the 
trigger for ULIRG activity in the local Universe. At low- 
z, ULIRGs are formed via the merger of two or more gas 
rich spiral galaxies. At high redshift however the galaxy for- 
mation processes were likely more diverse, with hierarchical 
buildup from many small dwarf galaxies or monolithic col- 
lapse of a large disk of gas forming a primeval galaxy both 
leading possibilities. 



7 CONCLUSIONS 

We have studied the properties of a sample of 41 local 
Ultraluminous Infrared Galaxies using archival optical 
and infrared photometry and advanced radiative transfer 
models for starbursts and AGN. Our conclusions are: 

1) All of the sample contain a luminous starburst, whereas 
about half contain a luminous AGN. The mean starburst 
fractional luminosity is 82%, and in ~ 90% of the sample 
the starburst produces more than half the total IR emis- 
sion. We conclude that the fraction of purely AGN powered 
ULIRGs in the local Universe is less than 2%. By combining 
our objects with other galaxy samples we find that starburst 
and AGN luminosities correlate over 6 orders of magnitude 
in total IR luminosity and over a wide range of galaxy types 
suggesting that a common physical factor, most plausibly 
the gas masses in the nuclear regions, govern both the star- 
burst and AGN luminosities. 

2) The starburst luminosity shows a strong positive corre- 
lation with total luminosity, as does the AGN luminosity, 
albeit less strongly. We however find no trend for increasing 
fractional AGN luminosity with increasing total luminosity, 
contrary to previous claims. We find that these claims were 
based on finding generally more luminous AGN in more lu- 
minous ULIRGs, rather than increasing AGN dominance 
with increasing total luminosity. 

3) We derive a mean starburst age range in ULIRGs of 
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1.0 x 10 7 — 4.0 x 10 7 years. Together with previous esti- 
mates for the lifetime of a ULIRG and an AGN, we find 
that most ULIRGs must undergo multiple starbursts dur- 
ing their lifetime. When combined with recent simulations 
of pair and multiple galaxy mergers we infer that mergers 
between more than two galaxies must be common in the 
ULIRG population. 

4) By comparison with previous results we find that the mid- 
IR F7.7/C7.7 line-continuum ratio gives no indication of the 
luminosity of the starburst in ULIRGs. We also find that 
Ft. 7/C7. 7 ratio gives no indication of the fractional AGN lu- 
minosity. We find therefore that a large value of F7.7/C7.7 
only indicates the presence of a moderately obscured star- 
burst, but gives no information on the presence or proper- 
ties of heavily obscured starbursts and AGN. As such, the 
-P7.7/C7.7 ratio on its own is not a reliable diagnostic of the 
power source in ULIRGs. 

5) The total 1.4GHz radio flux for the objects in our sample 
correlates strongly with the starburst luminosity but shows 
no correlation with the AGN luminosity. From this we infer 
that the radio-IR correlation in ULIRGs is due to star forma- 
tion, in line with previous results. Furthermore, we propose 
that the scatter in the correlation is due to a skewed IMF of 
the starburst and/or a relic relativistic electron population 
from a previous starburst event, rather than contamination 
from an obscured AGN. 

6) The rarity of luminous AGN and AGN dominated systems 
argues against a simple evolutionary model for ULIRGs in 
which they all evolve to become optical QSOs. Furthermore, 
although 'warm' ULIRGs generally contain more luminous 
AGN than do the 'cool' ULIRGs, there is no difference in 
fractional AGN luminosity between the 'warm' ULIRGs and 
the 'cool' ULIRGs. We find therefore that ULIRGs as a class 
do not evolve to become QSOs, but instead follow multi- 
ple evolutionary paths in transforming merging spirals into 
emerging ellipticals, and that only a few U LIRGs become 
optical QSOs, as suggested by iFarrah et alll20o"lh . 

7) By comparing our local sample to a sample of HLIRGs 
at z ~ 1 we find that AGN activity is much higher in the 
z ~ 1 sample. We infer that the two samples are distinct 
populations and postulate that different galaxy formation 
processes at high-z are responsible for this difference. 
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